Theoretical study of the conductance of ferromagnetic atomic-sized contacts 



M. Hafner, 1 J.K. Viljas, 1 ' 2 D. Frustaglia, 3 F. Pauly, 1 - 2 M. Dreher, 4 P. Nielaba, 4 and J.C. Cuevas 5 - 1 - 2 

1 Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 
2 Forschungszentrum Karlsruhe, Institut fur Nanotechnologie, D-76021 Karlsruhe, Germany 
3 NEST- CNR-INFM & Scuola Normale Superiore, 1-56126 Pisa, Italy 
J * Fachbereich Physik, Universitat Konstanz, D-78457 Konstanz, Germany 
5 Departamento de Fisica Teorica de la Materia Condensada, 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
(Dated: February 6, 2008) 

Recently, different experiments on the transport through atomic-sized contacts made of ferromag- 
netic materials have produced contradictory results. In particular, several groups have reported the 
observation of half-integer conductance quantization, which requires having full spin polarization 
and perfectly conducting channels. Motivated by these surprising results, we have studied theoreti- 
cally the conductance of ideal atomic contact geometries of the ferromagnetic 3d materials Fe, Co, 
and Ni using a realistic tight-binding model. Our analysis shows that in the absence of magnetic 
domains, the d bands of these transition metals play a key role in the electrical conduction. In the 
contact regime this fact has the following important consequences for the three materials: (i) there 
are partially open conduction channels and therefore conductance quantization is not expected, (ii) 
the conductance of the last plateau is typically above Go = 2e 2 /h, (iii) both spin species contribute 
to the transport and thus there is in general no full current polarization, and (iv) both the value 
of the conductance and the current polarization are very sensitive to the contact geometry and to 
disorder. In the tunneling regime we find that a strong current polarization can be achieved. 

PACS numbers: 73.63.Rt, 75.75.+a 



I. INTRODUCTION 



Metallic nanowires fabricated by means of scanning- 
tunneling microscope and break-junction techniques have 
turned out to be a unique playground to test basic con- 
cepts of electronic transport at the atomic scalei Usu- 
ally the conductance of these contacts is described by 
the Landauer formula G = Go^2 n T n , where the sum 
runs over all the available conduction channels, T n is the 
transmission for the nth channel, and Go = 2e 2 //i is the 
quantum of conductance. In the case of broken spin sym- 
metry one should include a sum over spin and replace Go 
by e 2 /h in the previous formula. It has been shown that 
the number of channels in a one-atom contact is mainly 
determined by the number of valence orbitals of the cen- 
tral atom, and the transmission of each channel is fixed 
by the local atomic environment^* 3 ^ Thus, for instance, 
a one-atom contact of a monovalent metal such as Au sus- 
tains a single channel, while for sp-like metals such as Al 
or Pb one finds three channels due to the contribution of 
the p orbitals. More importantly for the discussion in this 
work, in a transition metal such as Nb the contribution 
of the d orbitals leads to five partially open channels £££ 

In the last years a lot of attention has been devoted 
to the experimental analysis of contacts of magnetic 
matcriab . 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 
In particular, several groups have reported the obser- 
vation of half-integer conductance quantization in the 
last stages of the breaking of nanowires ,igiiLigii2iS2i2ii 
Moreover, some authors have shown that the conduc- 
tance histograms are very sensitive to a magnetic field. 
For instance, Ono et al. 11 reported that the peaks at 
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IGo and 2Go which are observed without magnetic field 
are joined by additional peaks near 0.5Go and 1.5Go for 
fields above 5 mT. These findings are rather surprising 
and challenge our present understanding of the trans- 
port properties of atomic-sized contacts. Half-integer 
conductance quantization would require to have fully 
spin-polarized contacts and conduction channels with 
perfect transparency. Neither of these properties are 
expected in the 3d ferromagnetic materials (Fe, Co, and 
Ni). In these transition metals both the majority spin 
and the minority spin electrons at the Fermi energy 
contribute to transport and, moreover, the partially 
occupied d orbitals are expected to give rise to partially 
open channels. Finally, it is impossible in practice for a 
magnetic field to lift the spin degeneracy in a metallic 
contact. A field of 10 T produces a Zeeman splitting on 
the order of 1 meV, which is extremely small in compar- 
ison with the bandwidth of even a single-atom contact, 
this being at least on the order of several electronvolts^ 
In this sense, in the absence of magnetic domains and 
magnetostriction, it is difficult to understand how an 
external magnetic field could significantly modify the 
conductance of an atomic contact. 

More recently, Untiedt et alM- measured the conduc- 
tance for contacts of several magnetic metals (Fe, Co, and 
Ni) using break junctions at low temperatures and under 
pyogenic vacuum conditions. They reported the absence 
of fractional conductance quantization, even when a high 
magnetic field was applied, which agrees better with 
the picture described above. Several recent model cal- 
culations support these findingsi 30 i 31 ' 32 i 33 i 34 i 35 i 36 i 37 i 38 i 39 
These studies have convincingly shown that neither con- 
ductance quantization nor full spin polarization are to be 
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expected in contacts of this kind. 

In spite of the coherent picture that is emerging, one 
still misses in the literature a comparative analysis of the 
ferromagnetic 3d materials (Fe, Co, and Ni) that clari- 
fies basic issues like which orbitals are relevant for the 
transport, the role of atomic disorder, or the dependence 
of the spin polarization of the current on the thickness 
of the contact. To fill this gap, we present in this work 
detailed calculations of the conductance of ferromagnetic 
atomic contacts of Fe, Co, and Ni. In our calculations 
we study ideal contact geometries using a realistic tight- 
binding model. We only consider situations where the 
contact region is a single magnetic domain. Our results 
indicate that for a few-atom contact of the three ma- 
terials one can draw the following general conclusions: 
(i) there is no conductance quantization, mainly due to 
the partially open conduction channels of the minority 
spin electrons, (ii) the last plateau has typically a con- 
ductance above Go = 2e 2 //i, (iii) the current is not fully 
spin-polarized and both spin species contribute to the 
transport, and (iv) both the value of the conductance and 
the current polarization are very sensitive to the contact 
geometry and disorder. The origin of all these findings 
can be traced back to the fact that the d bands of these 
transition metals play a very important role in the elec- 
trical conduction. This is in contrast with other physical 
situations such as tunnel junctions or bulk systems. Fi- 
nally, we find that in the tunnel regime, which is reached 
when the contacts are broken, the nature of the conduc- 
tion changes qualitatively and almost fully spin-polarized 
currents are indeed possible. 

The rest of the paper is organized as follows. In the 
following section we describe our tight-binding approach 
to compute the conductance of the ferromagnetic atomic 
contacts. Section HTT1 is devoted to the analysis of the re- 
sults of the conductance of representative one-atom thick 
contacts of Fe, Co, and Ni. Moreover, we include in 
this section a discussion of the conductance in the tun- 
nel regime. In Sec. lIVI we discuss the influence of atomic 
disorder on the conductance of single-atom contacts. Fi- 
nally, in Sec. we summarize and discuss the main re- 
sults. 



II. DESCRIPTION OF THE THEORETICAL 
MODEL 

Our goal is to compute the low-temperature linear con- 
ductance of atomic-sized contacts of the 3d ferromagnetic 
metals Fe, Co, and Ni. For this purpose, we use a tight- 
binding model based on the sophisticated parametriza- 
tion introduced in Ref. [ipj. Such tight-binding models 
have been successful in the description of electron trans- 
port in metallic atomic contactsi 2 ^^ Our approach fol- 
lows closely the one used in Refs. and we now 
proceed to describe it briefly. 

In our approach the electronic structure of the atomic 
contacts is described in terms of the following tight- 



binding Hamiltonian written in a non-orthogonal basis 

H= J] H icc,jpc\<x,<T t il3,a- (1) 

Here i,j run over the atomic sites, a, (3 denote the dif- 
ferent atomic orbitals, and a =t,| denotes the spin. 
Furthermore, H? a for i = j and a = (3 are the spin- 
dependent on-site energies, and for i ^ j the hopping el- 
ements, while H° a jg = for a ^ (3. In addition, we need 
the overlaps between the different orbitals, Siajp, which 
are spin- independent. We take all these parameters from 
the bulk parametrization of Ref. |40j |. which is known 
to accurately reproduce the band structure and total en- 
ergy of bulk ferromagnetic materials^ Notice that in our 
model there is no mixing of the two spin species, which 
means that, in particular, we do not consider spin-orbit 
interaction. The atomic basis is formed by 9 orbitals 
(3d, 4s,4p), which give rise to the main bands around 
the Fermi energy in Fe, Co, and Ni. It is important to 
emphasize that in this parametrization both the hopping 
elements and the overlaps are functions of the relative 
positions of the atoms, which allows us to study also ge- 
ometrical disorder. These functions have a cut-off radius 
that encloses atoms well beyond the 10th nearest neigh- 
bors in a bulk geometry for Fe, Co, and Ni. 

In order to compute the linear conductance, we apply 
a standard Green functions method 2 41 42 43 44] . For 
this we divide the system into three parts, the left (L) and 
right (R) leads, and the central cluster (C) containing 
the constriction. In this way, the retarded central cluster 
Green functions, Gq C , are given by 

G' CC {E) = [ES CC - H£ c - Hl(E) - E^)]" 1 , (2) 

Here Hcc anc ^ ^cc are the Hamiltonian and the overlap 
matrix of the central cluster, respectively, and are 
the self-energies, which contain the information of the 
electronic structure of the leads and their coupling to 
the central part of the contact. These self-energies can 
be expressed as 

Hi(E) = (Hcl - ES CL ) gl L (E) (Hl c - ES LC ) , (3) 

with a similar equation for E^(i?). Here, for example, 
H£, L is the hopping matrix connecting the central clus- 
ter C and the lead L, Scl is the corresponding block 
of the overlap matrix, and g^ L (E) is the retarded Green 
function of the uncoupled lead. Both infinite leads are de- 
scribed by ideal surfaces, the Green functions of which are 
calculated within the same tight-binding parametrization 
using the decimation technique described in Ref. |46j. 

In an atomic contact the local environment in the re- 
gion of the constriction is very different from that of the 
bulk material. In particular, this fact can lead to large 
deviations from the approximate local charge neutrality 
that typical metallic elements exhibit. We correct this 
problem by imposing charge neutrality on all the atoms of 



FIG. 1: (Color online) Bulk density of states (DOS) of Fe, Co, and Ni, resolved with respect to the individual contributions 
of 3d, 4s, and 4p orbitals, as indicated in the legend. Furthermore, the upper panels show the DOS for the majority spins and 
the lower ones the DOS for minority spins. The Fermi energy is set to zero and it is indicated by the vertical dashed line. 



the nanowire through a self-consistent variation of Hp C , 
following Ref. ^3 an d shifting both spin species equally. 

The linear conductance at low temperature can now 
be expressed in terms of the Landauer formula 

G=^T CT (i? F ), (4) 

a 

where T a (E) is the total transmission for spin a =|, j 
at energy E, and Ep is the Fermi energy. We also de- 
fine the spin-resolved conductances G a — (e 2 /h)T a (Ep), 
such that G = G| + G|. The transmissions are obtained 
as follows 

T a {E) = Tr[t CT (£)tt(£)] = ^T„ >CT (£), (5) 

n 

where t a (E) is the transmission matrix and T nt<T (E) are 
the individual transmission eigenvalues for each spin a. 
The transmission matrix can be calculated in terms of 
the Green functions G^c as follows 

t„{E) = 2 [Tl(E)} 1/2 Gq C (E) m(E)} 1/2 . (6) 

Here, T a L j R (E) are the scattering rate matrices given by 
n jR {E) = -Im[££ /R (£)]. 

III. CONDUCTANCE OF IDEAL 
SINGLE-ATOM CONTACTS OF FE, CO, AND NI 

The goal of this section is the analysis of the conduc- 
tance of ideal, yet plausible one-atom contact geometries 
for the three ferromagnetic metals (Fe, Co, and Ni) con- 
sidered in this work. In order to understand the results 
described below, it is instructive to first discuss the bulk 
density of states (DOS). The spin- and orbital-resolved 
bulk DOS of these materials around Ep, as calculated 



from our tight-binding model, is shown in Fig. ^ The 
common feature for the three ferromagnets is that the 
Fermi energy for the minority spins lies inside the d 
bands. This fact immediately suggests^ that the d or- 
bitals may play an important role in the transport. Oc- 
cupation of the s and p orbitals for both spins is around 
0.25 and 0.4 electrons, respectively. For the majority 
spins the Fermi energy lies close to the edge of the d 
band. The main difference between the materials is that 
for Fe there is still an important contribution of the d 
orbitals, while for Ni the Fermi level is in a region where 
the s and p bands become more important. The calcu- 
lated values of the magnetic moment per atom (in units 
of the Bohr magneton) of 2.15 for Fe, 1.3 for Co, and 
0.45 for Ni are reasonably close to literature values4£ 

We now proceed to analyze in detail the conductance of 
some ideal one-atom geometries, which are chosen to sim- 
ulate what happens in the last conductance plateau be- 
fore the breaking of the nanowires. First, we consider the 
one-atom contacts shown in the upper panels of Fig. |21 
These geometries are constructed starting with a single 
atom and choosing the nearest neighbors in the next lay- 
ers of the ideal lattice along the direction indicated with 
an arrow. In the case of Fig. [3 for Fe (bec lattice with 
a lattice constant of 2.86 A) the contact is grown along 
the [001] direction, for Co (hep lattice, lattice constant 
2.51 A) along the [001] direction (parallel to "c axis") 
and for Ni (fee lattice and lattice constant 3.52 A) along 
the [111] direction. The number of atoms in the cen- 
tral region has been chosen large enough, such that the 
transmission does not depend anymore on the number of 
layers included. Moreover, as explained in the previous 
section, the central region is coupled seamlessly to ideal 
surfaces grown along the same direction. 

Let us start describing the results for the Fe one-atom 
contact of Fig. EJa) . There we present the total trans- 
mission for majority spins and minority spins as a func- 
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FIG. 2: (Color online) Transmission as a function of energy for the three single-atom contacts of (a) Fe, (b) Co, and (c) Ni, 
which are shown in the upper panels. We present the total transmission (black solid line) for both majority spins and minority 
spins as well as the transmission of individual conduction channels that give the most important contribution at Fermi energy, 
which is indicated by a vertical dotted line. The blue, brown and violet dash-dotted lines of ri, t%, and T3 refer to twofold 
degenerate conduction channels. The legends in the upper graphs indicate, in which direction the contacts are grown. These 
contacts contain in the central region 59 atoms for Fe, 45 for Co, and 39 for Ni. The blue atoms represent a part of the atoms 
of the leads (semi-infinite surfaces) that are coupled to the central atoms in our model. 



tion of energy as well as the individual transmissions. 
We find for this particular geometry the spin-resolved 
conductances Gf = 3.70e 2 /h (f for majority spins) and 
G[ = 3.75e 2 /h (J. for minority spins), which results in 
a total conductance of 3.7Go- The conductance Gf for 
the majority spins is the result of up to 8 open chan- 
nels (with a transmission higher than 0.01), while for the 
minority spins there are 11 channels giving a significant 
contribution to Gj,. The large number of channels and 
consequently the high conductance, are partially due to 
the large apex angle of 71° of the pyramids. As a conse- 
quence of this, the layers next to the central atom couple 
to each other and give rise to a significant tunnel cur- 
rent that proceeds directly without traversing the central 
atom. On the other hand, the larger number of channels 
for the minority spins is due to the key contribution of 
the d orbitals that dominate the transport through this 
spin species, while for the majority spins the s and p 
orbitals are the more relevant ones. This fact, which is 
supported by the analysis of the local density of states 
(not shown here), is a simple consequence of the posi- 
tion of the Fermi energy and the magnitude of the spin 
splitting (see discussion of the bulk DOS above). 
We define the spin polarization P of the current as 



With this definition we find a value of P = —0.7% for 
the Fe one-atom contact of Fig. Efa) . In order to com- 
pare to the polarization of bulk we have calculated the 
transmission at the Fermi energy for a series of contact 
geometries where a bar of constant diameter bridges the 
two lead surfaces. When the diameter of the bar (or cen- 
tral region of the contact) is increased, the polarization 
grows continously and saturates at a value of P = +40% 
for a contact containing 219 atoms in 7 layers. This is in 
good agreement with the experimental value obtained us- 
ing normal-metal-superconductor point contacts^! No- 
tice that P can be quite different in an atomic contact 
as compared to bulk. This is because the conductance 
is not simply controlled by the DOS at the Fermi en- 
ergy, but the precise coupling between the orbitals in the 
constriction plays a crucial role. 

For the Co contact depicted in Fig. |2fb) the trans- 
mission is lower than for Fe, partly due to the smaller 
apex angle of the hep pyramids. In this case we find 
Gj = 1.57e 2 //i for majority spins and Gj_ = 1.21e 2 //i 
for minority spins, summing up to a total conductance 
of 1.4Go- There are 3 channels contributing to Gf and 8 
channels to Gf . As in the case of Fe, the larger number 
of channels for the minority spins is due to the position of 
the Fermi level and the resulting contribution of the d or- 
bitals for this spin. We also find that there is a small but 




FIG. 3: (Color online) The same as in Fig. |5|but for the geometries shown in the upper graphs, which contain a dimer in the 
central part of the contact. The two dimer atoms are at the bulk nearest neighbor distance from each other. 



non-negligible contribution of channels that proceed di- 
rectly without crossing the central atom. This explains, 
in particular, why one has 8 channels for the minority 
spins, although at most 6 bands (s and d) have a signifi- 
cant DOS at this energy. Turning to the current polariza- 
tion, we find a value of P = +13% for the Co one-atom 
contact. We also calculate the polarization for a series of 
Co bars with increasing diameter in hep [001] direction. 
As the diameter increases, the polarization decreases to 
a value of P = —41% for a contact containing five lay- 
ers of 37 atoms each, again in good agreement with the 
experiment. 51 Notice again that not only the magnitude 
of P for a one-atom contact can be quite different from 
bulk, but also its sign can be the opposite. 

Finally, the Ni contact shown in Fig. [2(c) exhibits 
conductances of Gf = 0.85e 2 //i for majority and Gf = 
1.80e 2 //i for minority spins, adding up to a total conduc- 
tance of 1.3Go- The Gf consists of 3 channels, due to 
the contribution of the s and p orbitals, and Gf contains 
6 channels, which originate from the contribution of the 
d orbitals. In this case we find a value for the polariza- 
tion of P = —34%. Once more we have investigated the 
polarization of bulk Ni in a series of large Ni bars in fee 
[111] direction. Interestingly, the polarization decreases 
from P = +3% for a contact of 28 atoms in four layers to 
P = —41% for a contact consisting of 244 atoms in four 
layers. 

Now we turn to the analysis of the geometries shown 
in the upper panels of Fig.[3J The difference with respect 
to the geometries of Fig. [JJ is the presence of a dimer 
in the central part of the contacts. This type of geome- 



try has frequently been observed in molecular dynamics 
simulations of atomic contacts of Al (Ref. [jsj) and Au 
(Ref. ^j|) and we also find them in our simulations of Ni 
contacts in the last stages of the breaking process^ 

Inserting a dimer in the geometries of Fig. [2J results in 
a larger separation of the pyramids to the left and right 
of the central atom and therefore in a weaker coupling 
between the layers next to the dimer. This is particu- 
larly important in the case of Fe. The resulting transmis- 
sion for the Fe contact with a central dimer is shown in 
Fig. [2(a), where one can see that only 3 channels remain 
for the majority spins, yielding Gf = 1.24e 2 //i, while 
for the minority spins 3 channels contribute to Gj = 
0.70e 2 /h. The total conductance is l.OGo and the polar- 
ization P = +28%. For Co the contact of Fig. [3(b) with a 
central dimer exhibits Gf = 0.90e 2 //i and Gf = 2.23e 2 //i, 
summing up to a total conductance of 1.6Go- The trans- 
mission is formed by 3 channels for the majority spins 
(with one clearly dominant) and 6 channels for minority 
spins and polarization is P = —42%. Finally for the Ni 
contact in Fig. [3(c) with a central dimer, a single channel 
contributes to Gf = 0.86e 2 //i and 4 channels add up to 
Gf = 2.66e 2 //i. This means that one has a total conduc- 
tance of 1.8Go, while the current polarization adopts a 
value of P = -51%. 

Beyond the precise numerical values detailed in the 
previous paragraph, we would like to stress the following 
conclusions from the analysis of Fig. [3| First, the trans- 
port contribution of the minority spins is dominated by 
the d orbitals, which give rise to several channels (from 3 
to 5 depending on the material). Second, for the major- 
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FIG. 4: (Color online) Conductance over tip separation D of similar geometries as in Fig. E| The conduct ance of majority spin 
(Gf, dashed lines and left scales) and minority spin (Gj, dash-dotted lines and left scales) is shown, together with the resulting 
current polarization (solid lines and right scales). 



ity spins there is a smaller number of channels ranging 
from 3 for Fe to 1 for Ni. This contribution is dominated 
by the d and s orbitals for Fe and only by the s orbitals 
for Co and Ni. The relative contribution and number of 
channels of the two spin species is a simple consequence of 
the position of the Fermi level and the magnitude of the 
spin splitting. In particular, notice that as we move from 
Fe to Ni, the Fermi energy lies more and more outside 
of the d band for the majority spins, which implies that 
the number of channels is reduced for this spin species. 
In particular, for Ni a single majority spin channel dom- 
inates. On the other hand, notice that the conductance 
values for the different contacts lie typically above Go, 
which is precisely what is observed experimentally^ 

So far, we have analyzed geometries for the so-called 
contact regime where the nanowires are formed. As 
shown above, in this case the contribution of the d bands 
makes it difficult to obtain large values of the current 
polarization. In this sense, one may wonder what hap- 
pens in the tunnel regime when the contact is broken. In 
order to address this issue, we have simulated the break- 
ing of the contacts by progressively separating the elec- 
trodes of the dimer geometries of Fig. [31 In this way, 
we have computed the conductance and the current po- 
larization as a function of the tip separation D and the 
results for the three materials are summarized in Fig. 0] 
With increasing D one enters the tunnel regime, which 
is characterized by an exponential decay of the conduc- 
tance. In the regime shown in the graphs, Fe does not 
yet exhibit an exponential decay. In contrast, the con- 
ductances for Co and Ni are well fitted by an exponential 
exp(-/3D) with (3 = 2.3 A" 1 and (3 = 1.9 A -1 , respec- 
tively. These values are in reasonable agreement with the 
WKB approximation^, which yields (3 = 2.2 A" 1 using 
a work function of 5 eV4£ Notice that deep in the tunnel 
regime for the three materials, the conductance for the 
majority spins largely overcomes the value of the minor- 
ity spin conductance. This results in positive values of 
the current polarization P and, in particular, for Co and 
Ni it reaches values very close to 100%. 

The origin of these huge values of current polarization 



in the tunnel regime is the following. In this regime the 
current is, roughly speaking, a convolution of the local 
densities of states on the tips weighted with the squared 
hoppings of the relevant orbitals of both electrodes. The 
hoppings between 3d orbitals decay faster with the sep- 
aration of the tips than the corresponding hoppings of 
the 4s orbitals. Due to this faster decay, the conduc- 
tion is then dominated by the s orbitals and since the 
on-site energy for the minority spins lies further away 
from the Fermi energy than the corresponding one for 
majority spin, the transmission through the latter one 
is much higher giving rise to a very large positive spin 
polarization P in the tunneling regime. 



IV. ROLE OF ATOMIC DISORDER 

In the previous section we have seen that the 3d or- 
bitals play an important role in transport. These orbitals 
are rather localized on the atoms and the energy bands 
that they give rise to have relatively flat dispersion re- 
lations. Therefore, one would expect the contribution 
of these orbitals to the transport to be very sensitive to 
the contact geometry. Indeed, in the previous section we 
have seen examples in which, by changing the structure 
of the central part of the contacts, one can even change 
the sign of the current polarization. Motivated by these 
results, we study in this section in a more systematic 
manner how disorder in the atomic positions influences 
the conductance of one-atom contacts. 

In order to simulate the role of disorder we have stud- 
ied the conductance of contacts in which the atomic posi- 
tions in the central cluster have been changed randomly 
using the geometries of Figs. [3] and [3] as starting points. 
In Fig. [S] we present an example of such a study, where we 
show histograms of the channel-resolved transmissions at 
the Fermi energy Tn^Ep) for both spins a constructed 
from around 3000 realizations of disorder for each con- 
tact. The amplitude of the random displacement in each 
direction was in this case 0.05 times the lattice constant. 
Similar results for contacts of the noble metal Au are 
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FIG. 5: (Color online) Histograms of transmission channels at Fermi energy, T n , a (EF), for 3000 perturbed realizations of ideal 
contact geometries of Fe, Co, and Ni. Panels (a)-(d) show histograms for contacts with a single central atom as in Fig. Hand 
panels (e)-(h) for contacts with a central dimer as in Fig. [21 In panels (a)-(c) and (e)-(g) results for ferromagnetic contacts are 
presented: the upper parts of the panels refer to majority spin channels and the lower parts to minority spin channels. Only 
channels that contribute more than 0.01 to transmission are displayed, and the histograms for smaller transmission values are 
in the front. The insets for the ferromagnetic materials show corresponding histograms for the current polarization P, where 
on the x axis P is given in %. Finally, panels (d) and (h) show comparison histograms for fcc-Au calculated with a similar set 
of geometries as for Ni. 



also shown for comparison. Moreover, for the ferromag- 
netic materials the insets show corresponding histograms 
of the spin polarization P of the current. 

Let us now discuss the main features of the transmis- 
sion histograms. First, they show that the number of 
channels obtained for the ideal geometries in the previ- 
ous section is robust with respect to disorder, although 



the transmission coefficients depend crucially on the pre- 
cise atomic positions. Second, for the minority spins one 
has a non- negligible contribution of at least 5 channels, 
which originate mainly from the d bands. For the ma- 
jority spins the number of channels is clearly smaller and 
is progressively reduced as we go from Fe to Ni. This is 
particularly obvious in the panels of the dimer structures, 
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where one can see that for Fe there are three sizable chan- 
nels and the contribution of the smallest two decreases for 
Co and Ni. As explained in the previous section, this is a 
consequence of the relative position of the Fermi energy 
in these three metals. For the latter case of Ni, one chan- 
nel clearly dominates the majority spin conductance, but 
second and third channels arc still present. Thus, unlike 
in the case of noble metals such as Au, which only have 
a single channel (see Fig. [5J, for ferromagnetic materi- 
als conductance quantization is not expected. Third, the 
peaks in the histograms for the ferromagnetic metals are 
much broader (especially for the minority spins) than for 
Au. This is due to the higher sensitivity of the d bands to 
the atomic positions, as compared to the s orbitals that 
dominate the transport in the case of Au. This higher 
sensitivity is a result of the anisotropic spatial depen- 
dence of the d orbitals. 

In addition, we have calculated the values of the cur- 
rent polarization P for each realization of disorder in the 
contacts. The resulting histograms can be found as in- 
sets in the panels of Fig. [S] The peaks in each histogram 
are centered around the polarization values of the corre- 
sponding ideal geometries of Sec. IIIII 

To end this section we would like to make the following 
comment. In this work we have analyzed the conductance 
of some ideal one-atom geometries and the influence of 
disorder. These types of calculations are very valuable 
to elucidate the nature of the electrical conduction in 
atomic wires. However, one has to be cautious in es- 
tablishing a direct comparison between such theoretical 
results and the experiments because of the lack of knowl- 
edge of the exact geometries realized experimentally. Ide- 
ally, the theory should aim at describing the conductance 
histograms, which contain the full experimental informa- 
tion without any selection of the data. This is precisely 
what we have done for Ni contacts in our recent work 
[50| and refer the reader to it for further details. 



ductance quantization, neither integer nor half-integer. 
Second, the current in these junctions is, in general, not 
fully spin-polarized. Third, the conductance of the last 
plateau is typically above Go- Finally, both the con- 
ductance and the spin polarization of the current are 
very sensitive to the contact geometry. The ensemble 
of these theoretical findings supports the recent obser- 
vations of Untiedt et al-r^ while it is in clear contra- 
diction with the observations of half-integer conductance 
quantization, 16 ! 17 ! 18 ! 19 ! 20 ! 21 

It is interesting to mention that in the tunnel regime, 
when the contacts are actually broken, the nature of the 
conduction changes radically. We have shown that in this 
case the transport is mainly dominated by the s orbitals 
and the spin polarization of the current can reach values 
close to +100%. 

We want to stress that in all our calculations we have 
assumed that the atomic contacts were formed by single 
magnetic domains. In this sense, it would be interest- 
ing to see how the conductance in these calculations is 
modified by the presence of domain walls in the contact 
region. The first theoretical studies^ii2iMi2& along these 
lines show that the presence of a domain wall cannot 
conclusively explain the appearance of huge magnetore- 
sistance values reported in the literature. 1 - 

Recently, the so-called anisotropic magnetoresis- 
tance has been observed in ferromagnetic atomic 
contactsii^iSL 2 ^ 2 ^ This effect, in other words the de- 
pendence of the resistance on the relative alignment of 
the current and the magnetization, stems from the spin- 
orbit coupling, and can give rise to a correction to the 
resistance on the order of 1% in bulk ferromagnetsS 
Although the correction can be bigger for atomic-sized 
contacts^ 8 ^ it is nevertheless expected to be a rela- 
tively small effect. The main ingredient that determines 
the conduction in the 3d ferromagnets is the electronic 
structure, which is what we have described in this work. 



V. CONCLUSIONS 

In this work we have presented a theoretical analy- 
sis of the conductance of one-atom thick contacts of the 
ferromagnetic 3d metals Fe, Co, and Ni. Our calcula- 
tions are based on a self-consistent tight-binding model 
that has previously been successful in describing the elec- 
trical conduction in non-magnetic atomic-sized contacts. 
Our results indicate that the d orbitals of these tran- 
sition metals play a fundamental role in the transport, 
especially for the minority-spin species. In the case of 
one-atom contacts, these orbitals combine to provide sev- 
eral partially open conduction channels, which has the 
following important consequences. First, there is no con- 
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